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Universal Properties of Two-Dimensional Boson Droplets 
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We consider a system of nonrelativistic bosons in two dimensions, interacting weakly via a 
short-range attractive potential. We show that for A'' large, but below some critical value, the 
properties of the A^-boson bound state are universal. In particular, the ratio of the binding energies 
of (A'^-l-l)- and TV-boson systems, Bjv+i/-B]v, approaches a finite limit, approximately 8.567, at large 
A'^. We also confirm previous results that the three-body system has exactly two bound states. We 



find for the ground state Bf ^ = 16. 522688(1)52 and for the excited state B^'^' = 1.2704091(1)5; 
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PACS numbers: 03.75.Nt, 05.10.Cc, 36.40.-c, 12.38.A's 

The recent experimental progress with ultracold 
atomic gases has revived the interest in weakly coupled 
quantum liquids. The ability to control the parameters of 
the systems make trapped atomic gases ideal laboratories 
where theoretical ideas can be checked versus experiment. 
One of the fundamental parameters that can be varied in 
experiments is the dimensionality of space. Both one- 
and two-dimensional Bose-Einstein condensates (BEC's) 
of sodium atoms have been studied in atom traps [1] . A 
one-dimensional condensate of ^Li atoms immersed in a 
Fermi sea of ^Li atoms was observed in [2]. In Ref. [3], 
a two-dimensional BEG of cesium atoms was realized in 
a gravito-optical surface trap. A two-dimensional boson 
system has also been realized in hydrogen adsorbed on a 
helium surface [4]. 

In this paper, we revisit the problem of weakly inter- 
acting bosons in two spatial dimensions (2D). While most 
previous theoretical studies were concerned with a Bose 
gas with repulsive interactions [5, 6], we focus on at- 
tractive interactions. In particular, we consider a self- 
bound droplet of iV(^ 1) bosons interacting weakly via 
an attractive, short-ranged pair potential. Our goal is 
to exhibit universal properties pertaining to large finite 
systems, which are not in the thermodynamic limit. 

We shall show that the system possesses surprising uni- 
versal properties. Namely, if one denotes the size of the 
A^-body droplet as Rn, then at large N and in the limit 
of zero range of the interaction potential: 



Rn+i/Rn -O.Mll, 7V>1 



(1) 



The size of the bound state decreases exponentially with 
A'': adding a boson into an existing iV-boson droplet re- 
duces the size of the droplet by almost a factor of three. 
Correspondingly, the binding energy of iV bosons Bn in- 
creases exponentially with N: 



Bn+i/B 



N 



.567, 



iV> 1 



(2) 



This implies that the energy required to remove one par- 
ticle from a A'^-body bound state (the analog of the nu- 
cleon separation energy for nuclei) is about 88% of the to- 
tal binding energy. This is in contrast to most other phys- 
ical systems, where separating one particle costs much 



less energy than the total binding energy, provided the 
number of particles in the bound state is large. 

To derive results independent of the details of the 
short-distance dynamics such as the ones quoted above, 
the A^-body bound states need to be sufficiently shallow 
and hence have a size Rn large compared to all other 
length scales in the problem. A similar reasoning has 
been used in 3D with much success [7]. The breakdown 
of universality is determined by the next largest length 
scale in the problem, the natural low-energy length scale 
£. Depending on the physical system, £ can be the van 
der Waals length lydw, the range of the potential or 
some other scale. For realistic systems, Eqs. (1, 2) are 
valid for large A^, but below a critical value. 



K A < A^c: 



0.931 ln(i?2/^) + . 



(3) 



At A^ = Acrit the size of the droplet is comparable to 
£ and universality is lost. If there is a large separation 
between i?2 and £, then A'crit is much larger than one and 
the condition (3) can be satisfied. 

Asymptotic freedom. — Our analysis relies strongly on 
the property of asymptotic freedom of 2D bosons with 
attractive interaction, so we will briefly review this prop- 
erty here. The system under consideration is described 
by the nonrelativistic Hamiltonian 



H 



(4) 



The bosons interact via an attractive, short-ranged pair 
potential —g5'^{x), with <? > 0. This choice can be made 
because at low energies, the true potential can not be dis- 
tinguished from a (5-function potential. For convenience, 
we will use the unit system h = m = 1; the factors of h 
and TO can be restored from dimensional analysis. In this 
unit system, g is dimensionless, and we assume g <C 1. 

In 2D, any attractive potential has at least one bound 
state. For the potential —gd'^{x) with small g, there is ex- 
actly one bound state with an exponentially small bind- 
ing energy. 



B2 - A^ exp (-47r/g) , 



(5) 
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where A is the uhraviolet momentum ciitofF (which is the 
inverse of the range of the potential). Equation (5) can 
be obtained by directly solving the Schrodinger equation. 
However, this method is not practical for a system with 
more than a few particles. 

Asymptotic freedom provides an alternative way to un- 
derstand Eq. (5). In 2D nonrclativistic theory, the four- 
boson interaction g{ip'^tfj)'^ is marginal. The coupling runs 
logarithmically with the length scale R, and the running 
can be found by performing the standard rcnormalization 
group (RG) procedure. The RG equation reads 



dg{R) ^ gHR) 
d\nR 2tt 



(6) 



Depending on the sign of g, we find two different behav- 
iors. For repulsive interactions with g < 0, the coupling 
becomes weaker in the infrared. For g > 0, the coupling 
grows in the infrared, in a manner similar to the QCD 
coupling [8, 9]. The dependence of the coupling on the 
length scale R is given by 



9{R) = 



i _ ± ln(A2i?2) 
5 47r ' 



(7) 



so the coupling becomes large when R is comparable to 



the size of the two body bound state 



-1/2 



This is in 



essence the phenomenon of dimensional transmutation; 
a dynamical scale is generated by the coupling constant 
and the cutoff scale. 

It is natural, then, that B2 is the only physical energy 
scale in the problem: the binding energy of three-particle, 
four-particle, etc. bound states arc proportional to B2. 
However, the iV-particle binding energy can be very 
different from B2 if N is parametrically large. We shall 
now argue that Bn increases exponentially with N. 

Stabilizing the size of the droplet. — We first try to use 
the variational method to estimate the size of the bound 
state. For a cluster of a large number of bosons, one can 
expect that classical field theory is applicable. We thus 
have to minimize the energy (4) with respect to all field 
configurations V'(x) satisfying the constraint 



(8) 



It is instructive to first minimize the energy with respect 
to the size of the bound state, and afterwards over all 
shapes of the wave function. We use the following trial 
wave function: 



V'(x) = 



N 



RV2^ 



(9) 



where r = |x| and f{r/R) is a function that describes the 

shape of the wave function. We assume that f{r/R) is 
nonzero when r/R ^ 1, but becomes small when r/R^ 



1. Thus R is the size of the droplet. To satisfy the 
particle number constraint (8), we should take 



C = Jdppfip). 



(10) 



The total energy is then obtained from Eq. (4) as 
AN B gN^ 

where A and B depend on the shape of the wave function, 
A = jdpp[f'{p)\\ B = jdppf\p). (12) 

As one can see from Eq. (11), in 2D both the kinetic 

and potential energies scale as R~^ . This seems to pro- 
hibit a stable bound state: for A^ < 2itAC /{Bg) the sys- 
tem expands to infinite size, while in the opposite regime, 
N > 2TTAC/{Bg), it shrinks to zero size. 

However, the above estimate is too crude because it 
fails to account for the logarithmic running of the cou- 
pling g. We will therefore replace g in Eq. (11) by the 
coupling at the length scale R: g ^ 9{R)- This pro- 
cedure goes beyond the naive mean-field treatment and 
captures all leading logarithms (by using the RG), but 
it does not take into account all 1/A^ corrections. We 
shall see that this is sufficient for finding the parametric 
dependence of Bn on A''. 

Once g is replaced by g{R) in Eq. (11), the energy 
E(R) has a minimum at a finite R. Indeed, in the limit 
i? — > the coupling becomes weak, g 0, and E{R) is 
dominated by the kinetic energy, which tries to make the 
system larger. In the opposite limit i? — > 00 the coupling 
constant becomes strong, and E{R) is dominated by the 
negative potential energy, which favors smaller R. 

To find the optimal R, we differentiate the energy (4), 
with g replaced by g{R), with respect to R. We find 



AC 



Ng{R)^ , Ng^iR) 
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B^O, 



where we have used Eq. (6). The solution is 



(13) 



(14) 



The coupling is 0{N~^), which implies the weak- 
coupling regime at large N. The 0{N~'^) correction is 
beyond the scope of the classical approximation. Using 
the RG running of the coupling constant (7), we find the 
optimal size of the droplet. 



RN = CRR2e^^[-^N 



(15) 



where Cr is some numerical constant of order 1, which 

cannot be found at the current level of approximation 
due to the 0{N^'^) uncertainty in Eq. (14). However, 
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FIG. 1: Numerical solution of Eq. (17) for the boundary con- 
dition /'(O) = /(oo) = obtained using the shooting method. 



we can already see that the size of the droplet decreases 
exponentially as a function of the number of particles. 

From Eq. (14) we see that the kinetic and potential 
energies cancel each other to leading order in 1/iV. For 
this reason, we can only estimate the energy to be 
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Ce Ce 



N 



(j2 ^2 exp 



2B 

AC' 



(16) 



(barring the possibility that there is a cancellation in the 
next-to-leading order in 1/N), but cannot compute the 
overall constant Ce- 

The shape of the droplet. — We now minimize the en- 
ergy with respect to the shape of the wave function 
f(r/R). Due to the exponential behavior of the energy 
as a function of N , the optimal shape is the one which 
maximizes the ratio B/{AC), where A, B, and C are de- 
fined in Eqs. (10, 12). This ratio is truly characteristic 
of the shape of the wave function — it is unchanged un- 
der the rescaling f{p) Ai/(A2p). The optimal shape 
of the wave function is therefore ambiguous up to this 
trivial rescaling. 

Taking the variation of B/{AC) over the /(p), we find 
that f{p) satisfies the equation 



fip) 
P 



where we have performed the rescaling 




(17) 



(18) 



The boundary condition on f{p) is /'(O) = /(oo) = 0. 
The solution can be found numerically by using, e.g., 
the shooting method. The solution, shown in Fig. 1, 
has a characteristic bell shape with /(O) ~ 2.206. For 
the shape given by the function f{p) solving Eq. (17), 
A=\B = C ^ 1.862, therefore B/IaC) « 1.074. Equa- 
tion (16) now can be written as 



B 



N 



C1B2C 



N-2 



(19) 





FIG. 2: The integral equation for the three-body amplitude. 
The single (double) line indicates the boson (full dimeron) 
propagators, respectively. 



where c ~ 8.567, but ci is still unknown. Equations (1) 
and (2) are also recovered. 

Equation (17) resembles the Hartree equation for the 
single-particle wave function. The above results can also 
be obtained in the Hartree approach, provided the run- 
ning coupling is used instead of the bare one. 

Three-body bound state. — We next describe our compu- 
tation of the binding energies of the three-body system, 
which can be calculated exactly. For this purpose, we 
use an effective field theory and work in the Lagrangian 
formalism. It is convenient to introduce an auxiliary field 
d = ilP' with the quantum numbers of two bosons (some- 
times called the "dimeron") [10, 11]. In terms of d and 
-0, the Lagrangian density corresponding to Eq. (4) reads 



2 



■d^d- 



(20) 



The boson propagator takes the usual nonrelativistic 
form i/{po — /2 + ie). It is not renormalizcd by inter- 
actions since all tadpole diagrams vanish in this theory. 
The bare dimeron propagator is simply a constant ~-2i/g. 
In the presence of interactions, it gets dressed by boson 
bubbles to all orders, leading to the full propagator: 



iA(po,p) = -i^ln 
9 



rpV4-po 



B2 



(21) 



where the bare coupling constant g will drop out of all 
observables in the end. The Feynman rule for the drptjj 
vertex coupling the dimeron to two bosons is ig. 

The three-body binding energies are determined by the 
homogeneous integral equation for the three-body bound 
state amplitude depicted in Fig. 2. The single (double) 
line indicates the boson (full dimeron) propagators, re- 
spectively, while the blob is the bound state amplitude. 
It depends on the total energy E and the relative momen- 
tum of the boson and the dimeron. The three-body bind- 
ing energies are given by those (negative) values of the to- 
tal energy E = —B3, for which the homogeneous integral 
equation shown in Fig. 2 has a nontrivial solution. The 
derivation of the integral equation using the Feynman 
rules given above proceeds as in the three-dimensional 
case [11]. There are only bound states if the dimeron 
and the third boson are in a relative S-wave. The for- 
mation of bound states in the higher partial waves is 
prevented by the angular momentum barrier. Projecting 



4 



onto the S-wave, we obtain the equation for the bound 
state amphtudc F{p): 



F{p) = f 
Jo 



AqdqF{q) ln[{3qyA + B3) / B2]' 



p2q2 



(22) 



The three-body binding energies can be obtained numer- 
ically to high precision by discretizing Eq. (22). 

We find exactly two three-body bound states: the 
ground state with binding energy Bf^ = 16.522688(1)^2 
and one excited state with S^^^ = 1.2704091(1)^2. (The 
numbers in parentheses give the error in the last digit.) 
The thrcc-body binding energies for a zero-range po- 
tential in 2D have previously been calculated by Bruch 
and Tjon [12] and Nielsen et al. [13]. Our results are 
consistent with the previous calculations but more pre- 
cise. Platter et al. have recently calculated the four- 
body binding energies for a zero-range potential in 2D 
[14]. They found exactly two bound states: the ground 

state with ^4"' = 197.3(l)i?2 and one excited state with 



B 



(1) 



25.5(1)52. 



The results for the ground state energies Bf^ and B, 
are what should be compared with the asymptotic for- 
mula (2). The ratio B^^ / B^ ~ 16.5 is almost twice 
as large as the asymptotic value (2), while the ratio 
B^^ /B^^ « 11.9 is considerably closer. Such deviations 
are expected for the small values of N we are dealing 
with. Note, however, that the ratio of the root mean 
square radii of the two- and three-body wave functions is 
0.306 [13], close to the asymptotic value (1). 

Conclusion. — We have evaluated parametrically the 
bound state energy of of N weakly attracting nonrela- 
tivistic bosons. Our results are obtained by minimizing 
the mean-field energy functional with a scale-dependent 
coupling. While our approximation is good enough to es- 
tablish the exponential behavior of the binding energy on 
the number of particles, it is not sufficiently accurate for 
evaluating the overall coefficient in front of the exponent. 
It would be valuable to develop a technique capable of 
doing so. 

We also have computed the binding energy of a system 
of three bosons. We have confirmed the previous finding 
that there are two bound states [12, 13] and improved the 
precision of the universal values for their binding ener- 
gies. One also would like to directly compute the ground 
state energy of the TV-body system for > 4 and com- 
pare the ratio B^^ / B^^_-^ with the asymptotic value. It 
would be also interesting to know the number of excited 
states of the A^-body bound system at large N. 

The many-particle bound state studied here is analo- 
gous in to the nontopological soliton, or Q-ball [15, 16] in 
relativistic quantum field theory. Nonrelativistic bosons 
in 2D provide an interesting example where the size of 
the nontopological soliton is stabilized by a pure quan- 
tum effect (the running of the coupling). 



?(0) 



An important question to explore is whether the result 
of this work can be extended to three-dimensional boson 
systems with large scattering length [17]. In contrast to 
the 2D case, this system displays the Efimov effect [18], 
and the three-body bound state energy B3 is an inde- 
pendent parameter. On the other hand, the four-body 
bound state energy can be expressed via B2 and B^ [19]. 
One would like to know if the binding energy of N bosons 
can be expressed in terms of B2 and Bs alone and find 
the large N behavior of Bn and the wave function. 

Finally, one should investigate the realizablity of self- 
bound 2D boson systems with weak interactions in ex- 
periments. According to the analysis of Ref. [14], the 
1/N corrections to Eqs. (1, 2) are small for A'' > 6. Using 
(3), this requires i?2/^ ^ 600. We are not aware of any 
physical system that satisfies this constraint. However, 
such a system could possibly be realized close to a Fesh- 
bach resonance where R2 can be made arbitrarily large. 
A interesting theoretical question is the dynamics of the 
droplet formation [20]. 
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